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We measure the dipion mass spectrum in X(3872) ^ J/i/) tt+tt” decays using 360 pb“^ of pp 
collisions at y/s = 1.96 TeV collected with the CDF II detector. The spectrum is ht with predictions 
for odd C-parity ^Pi, and ^Dj) charmonia decaying to J//: , as well as even C-parity 

states in which the pions are from decay. The latter case also encompasses exotic interpretations, 
such as a D°D*° molecule. Only the ^Si and J/ip p hypotheses are compatible with our data. Since 
^Si is untenable on other grounds, decay via J/ip pis favored, which implies C = +1 for the X(3872). 
Models for different J/ip-p angular momenta L are considered. Flexibility in the models, especially 
the introduction of p-u) interference, enable good descriptions of our data for both L = Q and 1. 

PACS numbers: 14.40.Gx 13.25.Gv 12.39.Mk 


The charmonium-like X(387^ stands as a major spec¬ 
troscopic puzzle. Its mass Q H, Q and what is 
known of its decays make assignments to the normal 
spectrum of cc states problematic @ . Its remarkable 
proximity to the D^D*^ mass—indistinguishable within 
uncertainties—has fueled speculations that it is a loosely 


bound deuteron-like D°D*° “molecule”, i.e. a {uc)-{cu) 
system mil Although a molecule is prominent among 
exotic interpretations, others have been proffered mi- 
Non-qq mesons are allowed within QCD, but an unequiv¬ 
ocal example remains elusive. Even as a conventional me¬ 
son, the A(3872) remains interesting as the cc spectrum 
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FIG. 1: Examples of “slices” in dipion mass m.^.^ of the J/ipn^n mass distributions, and their X(3872) fits. The raw yields 
Nx prior to efficiency corrections are quoted, and the arrow in the first panel marks the X(3872) mass. 


above the tfj{3770) is not well known. 

Insight into the X(3872) is offered by the dipion mass 
spectrum in X —> J/ip 7r''"7r“. Belle observed a preference 
for high TTTr-masses, contrary to expectations for triplet- 
D cc-states 0— the naive interpretation. Belle noted 
that X J14’ decay—isospin violating for charm- 

onium—produces high masses, and thus may be a hint 
for a D^D*^ molecule. Dipion spectra have been pub¬ 
lished and a preliminary analysis partially based 

on TTTT masses argues for a assignment [^, 

consistent with that expected for a D^D*^ molecule. 

At the Tevatron, large Ai(3872) samples are available, 
albeit with high backgrounds. Previously, we have mea¬ 
sured the X mass and confirmed the propensity for high 
TTTT masses . We also have made a preliminary measure¬ 
ment of the inclusive production fraction arising from b 
hadrons [m Here we measure the tttt mass spectrum. 

We use a sample of pp collisions at =1.96 TeV col¬ 
lected with the Collider Detector at Fermilab (CDF II) 
between February 2002 and August 2004. The detec¬ 
tor is described in detail elsewhere 113, and only the 
most relevant components for this analysis are summa¬ 
rized here. The central tracking system is immersed in 
a 1.4 T solenoidal magnetic field for the measurement 
of charged particle momenta pt transverse to the beam¬ 
line. It is composed of six layers of silicon-strip detectors 
(LOO and SVX ^3) surrounded by an open-cell drift 
chamber called the Central Outer Tracker (COT) M 
The active volume of the COT is a 3.1 m long cylinder 
with 8 superlayers of 12 wires each. The outermost de¬ 
tection system is planes of multi-layer drift chambers for 
detecting muons [l^. The Central Muon system (CMU) 
covers It^I < 0.6, where pseudorapidity p = — ln[tan(0/2)] 
and 9 is the angle of the particle with respect to the di¬ 
rection of the proton beam. Additional chambers (CMX) 
extend the muon coverage to \r]\ = 1.0. 

A dimuon trigger is used to obtain a. J/4 ^ 
sample. At Level 1 of a three-level trigger system, the eX- 
tremely Fast Tracker (XFT) Ea uses COT information 


to select tracks based on pT- XFT tracks with pr > 1.5 
(2.0) GeV/c are extrapolated to the CMU (CMX) cham¬ 
bers and compared with the positions of muon-chamber 
tracks. The event passes Level 1 if two or more XFT 
tracks are matched to muon tracks. Opposite-charge and 
opening-angle cuts are imposed at Level 2. At Level 3, 
full COT tracking information is used to reconstruct 
candidates. Events with candidates from 2.7 to 
4.0 GeV/c^ in mass are recorded for further analysis. 

This analysis 0 is based on an integrated luminosity 
of 360 pb“^. Candidate selection follows Ref. 0 with 
two exceptions (see below). After constraining p~ 
candidates to a common vertex, the dimuon mass must 
be within 60 MeV/c^ (~ 4 standard deviations) of the 
J/4 mass ^ 3 . This one degree-of-freedom (DoF) fit 
must have < 15. Pairs of charged tracks, each with 
Pt > 0.4 GeV/c and assumed to be pions, are fit with the 
p^p~ tracks to a common vertex. In this fit, the dimuon 
mass is constrained to the J/4 mass, and we demand 
X^ < 25 (6 DoF). We reduce combinatorial backgrounds 
by requiring pt{J/4) > 4 GeV/c and Ai? < 0.7 for both 
pions, where Ai? = -^/qA^p+^A?))^, Ac/) is the differ¬ 
ence in azimuthal angle between the J/'c/’tttt system and 
the pion, and Ary is the difference in pseudorapidity. The 
mass range for the sample includes both A(3872) and 
4i‘^S) signals j3|; the latter is used as a control sample. 

We depart from Ref. by dropping a cut on the num¬ 
ber of candidates allowed per event. This removes a pos¬ 
sible bias and improves the X signal at high tttt masses. 
We also add fiducial criteria: pt{J/4 > 6GeV/c 
and |r 7 (J/r/)7r7r)| < 0.6. This eliminates the region of 
rapidly changing efficiency and sacrifices ~ 25% of the 
4i‘2S) yield, leaving 11500 ± 220'!/)(2S') mesons. We have 
1260±130A(3872) candidates for > SOOMeV/c^ @]. 

To extract dA/dmjr-n- spectra, we divide the sample 
into “slices” of and fit each J/r/;7r+7r“ mass distri¬ 
bution for the signal per slice. The J/r/irr+Tr” mass fits 
use a Gaussian for the A(3872) signal and an exponen¬ 
tial times power law for the background. We also fit the 
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'ip(2S) control signal in the same way, but use two Gaus- 
sians for the better defined shape arising from the larger 
'0(25') signal. As slices may have small signals—or none 
at all—we inhibit the fit from latching onto fluctuations 
by fixing the position and width of the signal to values 
from full-sample fits. Sample slices are shown in Fig.^ 

The fitted dN/dniTrTr yield is corrected for detector and 
selection efficiencies determined by Monte Carlo simula¬ 
tion. Only the efficiencies relative to other m,r-7r slices 
are needed. An important input affecting the efficiency 
is the production pT spectrum. For the Ai(3872), there is 
no a priori model, and we rely upon data. The generated 
spectra, exponentials raised to a quadratic polynomial in 
Pt, are adjusted until the simulation, after detector and 
reconstruction effects, reproduces the respective Ar(3872) 
or 0(25) spectra of the data. We use a parameterization 
of the well-known m,rir shape of the 0(25) for both states. 

Uncertainties on our spectra are dominated by statis¬ 
tics, but we examined two sources of systematic effects: 
the fits for signal yields and the efficiency corrections. 

To check the yield stability, we changed the width of 
the J/07r+7r“ fit range of 200 MeV/c^ by ±50 MeV/c^, 
altered the signal and background models, and allowed 
the signal mean and width to float. We saw no bias in 
the yields, but nevertheless allotted an uncertainty based 
on the statistical precision to which we could observe 
one: 3.6 0(25) and 8.4 ^(3872) candidates per slice. 
The highest three X-slices—highest two for the 0(25)— 
are treated specially for effects near the upper kinematic 
limit: the background begins to turn-on under the signal, 
and resolution effects can distort the signal shape in the 
mass fit. Yield systematics are assigned to these slices 
based on variations in the fit model for these issues. 

The other type of uncertainty is from the efficiency cor¬ 
rections. The 0(25) model for dA^/dm^rTr is inexact. We 
assign an uncertainty based on phase space as an alter¬ 
nate shape—including re-tuning the pr spectrum. The 
ratio of the alternate correction to the nominal one quan¬ 
tifies the change in shape of the efficiency when switching 
from the 0(25)-like dN/dm-^T^ to phase space. The ratio 
of efficiencies for the 0(25) varies by ^ 3% over the main 
region of interest, > 360 MeV/c^, and by ~ 2.5% for 
the ^(3872) above 570 MeV/c^. For the uncertainty in 
the meson px spectra, we use alternate spectra one stan¬ 
dard deviation steeper and shallower in their px fall-off 
based on the errors from the px spectrum fit to the data. 
We again take the ratio of the new efficiencies relative to 
the nominal shape to quantify the uncertainty. For the 
mass ranges of interest the 0(25) variation is almost 3%, 
but less than 1% for the X. The px spectrum of the X 
is more poorly measured than for the 0(25); but with 
higher dipion masses, the X suffers smaller variations in 
efficiency, and thus a smaller uncertainty. 

The dN/dm-rnc spectrum for our ip{2S) control signal, 
after corrections and including systematic uncertainties, 
is shown in Fig. |21 with a scale preserving the raw fitted 
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FIG. 2: The CDF dipion mass spectrum for the 0(25) with 
a fit of a QCD multipole expansion calculation for ®5i 1^ . 
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FIG. 3: The dipion mass spectrum for the X(3872), and fits 
to various hypotheses (see text). The fitted curve for the 
model is scaled up by a factor of 5 for better visibility. 


yield of 11500 candidates. It agrees well with results 
from the BES Collaboration using a sample of 20 000 
events This is reflected by the mutual agreement 
in fits to a QCD multipole expansion model j2]J. BES 
obtained —0.336 ± 0.009 ± 0.019 for this model’s single 
shape parameter “B/A,” whereas our fit yields —0.342 ± 
0.022 (6.9% fit probability). The systematic uncertainties 
are incorporated in this and later Y-fits, including the 
TO^r-TT correlations in the efficiency uncertainties. 

The X{3872) dipion spectrum is shown in Fig.|3 We 
fit our data with multipole expansion calculations for 
C-odd cc options jl^. The ^Dj states are a natural 
choice for the X accor ding to potential models |^. A 
fit of ^Dj J/0 7r+7r“ [2l| is unacceptable with a of 
113 for 14 DoF. The ^Pi J/07r+7r“ fit is worse 
/DoF = 146/14). The 0(25) spectrum is similar to 
that of the X, and indeed, our ^5i fit to the X has 
a 28% probability. However, no new ^5i cc-state can be 
near 3872MeV/c^ as the 0(35) lies at ~4040 MeV/c^ |^. 

The above C-odd states produce dipions, to lowest L 
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between the pions, with = 0 ++. C-even states yield 
1 isovector dipions, which we associate with the 
Isospin conservation suppresses cc decays to J/'i/' P- Thus, 
this mode is seen as suggestive of a D^D*° molecule w 
Even as charmonium, however, the X may break isospin 
by coupling to D^D*° due to its close proximity in mass. 

We model X ^ Jj'^ p hy relativistic Breit-Wigner 
multiplied by phase space and generalized for a J/V'-p 
system of angular momentum L. That is, dN/drriT^T^ oc 
(A:*)|i?pp, where fc* is the J/'0 momentum in 
the X rest-frame. Bp oc p/[w?p — rri^^ — irUpTp\, 

rp(m,r7r) = To [g*/go]^ [wp/m^,r][/ip(gH.)//ip(go)]^, q* is 
the TT momentum in the p rest-frame, go = g*(wp), Wp 
is 775.8 MeV/c^, and Tq is 146.4 MeV. The fLiip) are 
foi{p) = 1 and fii{p) = (1 + R^iP^)~^ [H, where Ri is 
a radius of interaction for meson “i.” The Ri are poorly 
known. A common value for light mesons is 0.3 fm, but 
for D mesons larger values like 1 fm are often taken . 
We use these respective values for Rp and Rx- 

Fits with this p model are shown in Fig. O for L = 0 
and 1. Higher L softens the fall-off at the high kinematic 
limit, worsening the agreement: the fit probability goes 
from 55% for L = 0 down to 7.7% for L — 1 |^. The 
P-wave fit is somewhat disfavored, but the results are 
sensitive to Rx and Rp. The latter probability can be 
increased by lowering Rp or raising Rx- 

Another modeling uncertainty is the effect of p-uj inter¬ 
ference. Belle reports evidence for A ^ J/'i/i 7 r+ 7 r“ 7 r° and 
interprets it as decay via a virtual uj. As such, they mea¬ 
sure the ratio of Jhhuj to T/'f/’P branching ratios P 3/2 to 
be 1.0±0.4±0.3 |^. The rate of a; ^ tt+tt" is normally 
negligible, but its interference effects may not be. 

We generalize |PpP to |ApPp -|- in 

dN 2 TT/dmTTT^, where Ap and are (positive) decay am¬ 
plitudes via p and w, and (j) is their relative phase. Using 
dN^^/dm^^ cx for 7 r+ 7 r“ 7 r° P 3/2 de¬ 

termines |A,^/Ap| given cf). We take a ((> of 95°, the value if 
the only phase is from uj —> 7 r''" 7 r“ decaying via p-LU mix¬ 
ing 1^ . Similar phases are seen in e+e“ —> tt+tt" [3(i| . 
The UJ fraction is small (< 10%), but interference is con¬ 
structive and contributes ^23% for both L, preferentially 
at high masses. Fits with this model are shown in Fig. 2] 
along with the breakdown into interference and “pure” p 
and UJ parts. The probability is 19% for the S fit, and 
53% for the P fit. The results are not critically depen¬ 
dent on 'Rzj'i- probabilities remain above 7% over a ±1 
standard deviation span of Belle’s P 3/2 for both L val¬ 
ues. The P fit is sensitive to (j) and Rx, as is shown in 
the inset. We conclude that there is ample flexibility in 
models ot X ^ J/tjjp oi either L to describe our data. 

In summary, we measured the dipion mass spectrum 
in A(3872) ^ J/’ip . Our spectrum is inconsistent 
with calculations for ^Pi and ^Dj charmonia. A good fit 
is obtained for X —> J/ip p^, an interpretation supported 
by recent evidence for the C-even decay X J/ijJl 113 ■ 
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FIG. 4: A blow-up of the A(3872) spectrum with J/ipp 
fits which include p-u interference (95° phase) with relative 
amplitudes set by 77.3/2 = 1.0. Fits for both L — 0 (lines) and 
1 (shaded regions) are shown, along with their decomposition 
into p, u, and interference terms. The inset shows L = 1 fit 
probabilities as a function of (f> and Rx in 5% contours. 



Our data are compatible with both S- and P-wave J/'i/' p 
decays, where in the latter case, this is partly due to 
modeling uncertainties. The P ht benefits from con¬ 
structive p-LJ interference at levels implied by the rate 
of A ^ J/V: 7 r+ 7 r“ 7 r°. The T/'i/'P interpretation does not 
by itself distinguish between C-even charmonia (e.g. I'*"'*" 
or 2 ^) and exotic options like an 1++ D^D*^ molecule. 
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